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Abstract 

The infinite divisibility of the squared Gaussian Process with 
non-zero mean is dependent on the signs of the row sums of the 
inverse of its covariance matrix. We give an explicit formula de- 
scribing this relationship. Let the Gaussian vector have a critical 
point «o if its squares are infinitely divisible for all a < ao and are 
not infinitely divisible for all a > ao- We give the upper bound 
for the critical point of a non-zero mean Gaussian vector. 

1 Introduction 

Let G = (r)x,r)2, ■■■■r) n ) be an n-dimensional Gaussian vector. We say that 
G has infinitely divisible squares (or G 2 is infinitely divisible) if for any 

m e N 

m 

G^= (itUV-.^ = 5> (1-1) 

1=1 

where Zi are independent identically distributed n- dimensional random 
vectors. 

Let N be an n x n matrix such that iVy = for all i j; N^i = ±1 
for alH. iV is called a signature matrix. 
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A non-singular matrix A is called an M-matrix if A^j > for all 
1 < i,j < n and Aij < for all i ^ j. Note that Aij must be 
positive. 

A symmetric matrix A is irreducible if it can not be written as a 
direct sum of square matrices. 

In what follows all covariance matrices considered are irreducible. 
The 1 denotes the n-dimensional vector with all entries equal to 1. The 
following Theorem is due to Griffiths and Bapat, see also [H Theorem 
13.2.1] and |4J Theorem 1.1]. Theorem 11.11 completely characterizes zero- 
mean Gaussian processes with infinitely divisible squares. 

Theorem 1.1. Let G = {G^G^-, G n ) be a mean zero Gaussian ran- 
dom variable with a strictly positive-definite covariance matrixT = = 
{EiGiGj)}. G 2 is infinitely divisible if and only if NT~ 1 N is an M- 
matrix for some signature matrix N . 

We are interested in Gaussian vectors with infinitely divisible squares. 
We can therefore restrict our attention to Gaussian vectors with positive 
covariance. 

The necessary and sufficient conditions for infinite divisibility of Gaus- 
sian squares with uniform non-zero mean - 

(G + al) 2 := ((£?! + a) 2 , (G 2 + a) 2 , ...) 

was posed and answered by Eisenbaum and Kaspi, [5J. Their result was 
developed to include the non-uniform mean Gaussian squares by Marcus 
and Rosen. Let c = (ci,C2-..) be any real n-dimensional vector. The 
question about infinite divisibility of (G + etc) 2 := ((G\ + ctci) 2 , (G2 + 
ac-i) 2 ...) for all a G R is answered by Theorem 1.3, [1]: 

Theorem 1.2. Let G be a mean zero Gaussian vector with strictly pos- 
itive definite covariance matrix T, N(0, b) denote the Gaussian random 
variable with mean and variance b. Let c = (ci, Oz---) and C be an n x n 
matrix with = Ci, Ci j = for i ^ j. The following are equivalent: 

1. (G + ca) 2 is infinitely divisible for all a 6 R. 



2 



2. There exists some b G R, b 7^ such that for rj = N(0,b), inde- 
pendent of G, (G\ + Ci?7, G 2 + c 2 r], ...G n + c n r],r]) has infinitely 
divisible squares. In this case (Gi + ci?y, G2 + C2^, ...G n + c n r],r]) has 
infinitely divisible squares for all b G R. 

3. CT _1 C is an M-matrix with non-negative row sums. 

Corollary 1.1. Let G, T, C be as above. When 1, 2 or 3 of Theorem 
hold, 

1- Ci t i 7^ 0, for all i. 

2. Fij^O, 

The natural question that could be asked at this point (and which 
was posed by Marcus and Rosen, see [2] and [5]) is about the existence 
and characterization of the vectors G for which (G + ac) 2 is infinitely 
divisible for some but NOT ALL a. We shall concentrate on the case 
c = 1. To put this question formally we need 

Definition 1.1. Let G be a Gaussian vector. Let «o be a real number 
< a < 00. We say that a is a critical point of G if 

• {G + al) 2 is infinitely divisible for all \a\ < a^, 

and (G + al) 2 is not infinitely divisible for any \a\ > cto 

• when (G + al) 2 is infinitely divisible for all a then G has no critical 
point. 

Marcus and Rosen have shown in [5] the existence of a critical point 
< ao < 00, for all two-dimensional Gaussian vectors with infinitely 
divisible squares. The question about existence of a critical points for 
the vectors of higher dimension remains open. 

We explore the connection between the properties of the inverse co- 
variance matrix and the infinite divisibility of (G + a) 2 for different 
a. We obtain some upper bounds on a critical point for Gaussian n- 
dimensional vectors G with n > 2 and also the proof of Theorem 11.21 
( part 3 1) which does not invoke the Markov process theory (see 
Section 4). In particular, the original proof of II. 21 uses the Second Ray- 
Knight Theorem to establish an isomorphism between (G + c) 2 and the 
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sum of two independent infinitely divisible random vectors. We obtain 
an elementary proof of the same result using the well known result due 
to Griffith and Bapat, (see [3]), Theorem 11.31 and Lemma [1.21 

We employ the necessary and sufficient conditions for infinite divis- 
ibility of a random vector due to Feller; (see [2J, Ch. XIII. 4), [I], [21 
13.2.2], and the application therein to Gaussian vectors. We follow the 
presentation of Marcus and Rosen in [3] Lemmas 13.3.1 and 5.2.1], and 

Lemma 1.1. J3| Lemma 13.2.2] Let G be a Gaussian vector with positive 
definite covariance T, a G R. For any vector A G R n denote by \I/g(A) 
the Laplace transform of (G + a) 2 and let S be a diagonal n x n matrix 
with entries G (0, 1]. Let t > 0. Let A be a diagonal matrix such that 
Aj j = Aj for all i. For any A G R n we can write 

A = t{I-S) (1.2) 

for any t sufficiently large and some S , where S is a diagonal n x n 
matrix with Sa G [0, 1]. Let 

= log tf G (A) 

and suppose that $(£, S) has a power series expansion about S = 0. Then 
(G + a) 2 is infinitely divisible if and only if for all t sufficiently large all 
coefficients of this expansion are non-negative, except the constant term. 

We obtain the series expansion for $(i, S) for general a in a form that 
makes it relatively easy to analyze the sign of the coefficients. 
Let 

Q(t) = [I + (tT)- 1 ]- 1 = /-(/ + tT)- 1 . (1.3) 
Let D G R n such that 

n 

fc=l 

i.e Di denotes the i-th row sum of r -1 . 
Let 

where 

$ 1 (t,S) = $(t,S)\ a ^ 
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The following series expansion for $i(t, S) is the result of Bapat and 
Griffith: 



» 1 (t,fl = l(iog|j-q| + £ tra "< qs > m ). (1.5) 



?Ti=l 

We extend their result to include a e R: 
Theorem 1.3. Under the notation as above, 

and for all a G R 



$ 2 (so = \ ( -i-r- 1 -! 7 + d - ( I £(gsr ) ■ d t ) (i + o(i/t)). 

(1.6) 



771=1 



Since (G + a) 2 is infinitely divisible if and only if the series expansion 
for $(5) has non- negative coefficients for all non-constant terms, it is 
clear from Theorem 11.31 that the signs of coefficients in the expansion of 
<& 2 (S)- and hence the infinite divisibility of (G + a) 2 will depend on the 
signs of the entries of Q and D. In particular, since the coefficients of $2 
depend linearly on a 2 , it is clear that if (G + a) 2 is not infinitely divisible 
for a particular a then it is not infinitely divisible for all /3 > a, and so the 
Definition 11.11 makes sense. The following Lemma establishes the positiv- 
ity of the entries of Q, which makes possible a relatively straightforward 
analysis of the coefficients of $(S). It also provides an alternative, ele- 
mentary proof of Theorem 11.21 which does not invoke the Isomorphism 
Theorems and Markov Processes. 

Lemma 1.2. For allt sufficiently large, 

Q>0 

Theorem 11.31 and Lemma 11.21 allow us to obtain an upper bound for 
the critical point of G: 
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Theorem 1.4. Let G be a mean zero Gaussian vector with infinitely 
divisible squares and covariance matrix T (hence T~ l is an M matrix). 
Let 

V = {{i,j) : A^<0}. (1.7) 
7/D ^ and «o is a critical point of G, then 



a < inf 

v 




An immediate consequence of Theorem 11.41 is 



Corollary 1.2. 

• If for some pair (i,j) e V Yj^ = 0, then a = 0. 

• ifV = ® then G has no critical point. 

Note that for the T> to be non-empty T -1 matrix has to have at least 
one negative row sum. Hence the second result of the Corollary is the 
statement ( 3 1) of Theorem 1 1.21 The first is the new result that allows 
us to construct examples of Gaussian vectors with infinitely divisible 
squares and zero critical point (see Section 7). We prove Theorem 11.31 
and Lemma [1.21 in Section 2. Theorem 11.41 and Corollary 11.21 are proved 
in Section 3. Section 4 contains an elementary proof of Theorem 11.21 
Various applications of Theorem II .31 are given in Sections 5 and 6. Section 
7 contains some numerical examples. 



2 Proofs of Theorem 11.31 and Lemma 11.2 

We use the following 
Lemma 2.1. 



$ 2 (S) 



aH 



1 • [(Q - I) + (I - Q- 1 ) (QS) m (Q - I)} ■ l T (2.1) 



m=l 
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Proof. By [3J, Lemma 5.2.1, 

v& G (A) = det(J + TA)- 1 / 2 • (exp(yl[Af A - A]1 T )). (2.2) 

where f = (T^ 1 + A) -1 . 
Let 

A = t(I - S) 

for sufficiently large t. 
Then 

$ 2 (t, 5) = y l[ALA - A]1 T =: y P(t, 5). (2.3) 

Let 

Q(t) = [i + (try 1 ]- 1 = /-(/ + try 1 . (2.4) 

For all sufficiently large t the expression in ( 12. 4ft renders itself to repre- 
sentation as an absolutely convergent geometric series: 

00 T-l — 1 

w = D- 1 ) w (— )" (2 - 5) 

u=0 

Hence for all 

<!•,-■ Qij ^M-- (2-6) 

We proceed to find the series expansion for $2(^5 S) about S = in terms 
of Q(t). Henceforth the parameter t will be suppressed in the expressions 
for Q(t) and S) for i = 1, 2. 

To obtain (12. ip consider (12. 3p and substituting for A using (11.21) write: 

Af a - a = t{i - ^[(tr)- 1 + / - sy\i -3)- t(I - S) (2.7) 

Hence 

p(s) = ti-{{i-s)[{try 1 + i-sy 1 {i-s) + {s-i)}-i T 
= a ■ {(i - s^q- 1 - sy\i -s) + (s- 1)} ■ i T 

= tl-{(I-S)[I-QS]- 1 Q(I-S) + (S-I)}-l T (2.8) 
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Since det(QS) < 1 for all t sufficiently large, 

oo 

[i-QS]- 1 = J2(QS) m - 

m=0 

Therefore 



(I-S)[I-QS]- 1 Q(I-S) + (S-I) (2.9) 

= (/-£)(£; (Qsr ) q(i -s) + (s-i) 



\m=0 



/ \m=0 / 

oo \ / oo 

E(W" +5" E(^) m ) +(^- J )- 



\m=l / \m=l 



For i — 0, 1 we write 

to see that the last line of ( 12. 9 p 

= (Q -/) + (/ - Q- 1 ) ( £ (QS) m ) (Q - J). 
This gives us (12. ip . 



\m=l 



□ 



Proof of Theorem 11.31 

Let 

n 

A^r-i; (2.10) 

fe=l 

i.e denotes the i-th row sum of V" 1 . 
By (123): 

/ - Q- 1 = I - (/ + (tr)- 1 ) = ZL2, (2.11) 
8 



and 

Q- 1=^(1 + 0(1/0). (2-12) 

Let D be a diagonal matrix such that D^i = Di. Now substituting this 
in (I2.1D we see that 



1 ■ [t(I - Q- 1 ) £ (QS) m (Q - I)} ■ 1 T (2.13) 

m=l 

oo 

= - 1 ■ r-^QS'JT-^l + 0(1/*)) • 1 T 

m=l 

_ oo n 

= i E E {r- 1 ^^- 1 }^^ + o(i/t)) 

m=l i,j=l 
1 oo n n 

= - t E E E r S(^ r S (i + o(i/t)) 

m=l i,j=l k,l=l 
1 oo n n n 

= ~ t E E (E r S)(W(E r Tj) (i + o(iA)) 

m=lfc,Z=l i=l i=l 
_ oo n 

= tEE(WA(i + o(iA)) 



m=l fc,/=l 



1 oo n 

m=l fc,/=l 
1 n / oo \ 

k,l=l \m=l / 

Now substitute this into (12. ip to get (II .61) . □ 
Proof of Lemma 11.21 

This is the direct consequence of the following 

Lemma 2.2. Let T~ l be an n- dimensional M matrix and assume that 
for some i j, T^j = 0. Let 

k = mm{l > 1 : rH ^ 0} 



Then 



k<n-l (2.14) 

and 

(r-% = (-i)1(r-%|. (2.15) 

Proof To show (12.151) make the following claim: if for all 1 < u < k, 
r~J = 0; then for any h < k, any term of the form 

V.^Vr^r- 1 ....r- 1 . 

/ j / j i,n r 1 ,r 2 ---- J - r h ,j 

equals 0. 

For k = 2; 

r=l 

since all summands are the products of two negative factors, hence pos- 
itive - the only potential negative summands are of the form T^jF^j, 
which is zero here, hence, if (r~ 2 )jj = 0, then for all 1 < r < n , 
ryirY- = - so the above claim holds. 

In general, for arbitrary k, suppose the claim holds for all I < k; 
Trj = 0. Then 

( r ~% = E E E r ^ F i 3 KLr (2.17) 
Suppose a particular summand contains one or more factors of the form 

r- 1 . 

Rearranging the summation order we get: 

VVVr 1 r- 1 ...r- 1 fV...Vr~' r~' ...r~ l ■". 

/ j / j / j r-L.ri 1 - r2,r 2 --- J - r u ,r u [ / 4 / „ x i,r u+1 J - r u+ i,r u +2 r k _± ,j 
ri T2 r u r u+ i r k _ 1 

(2.18) 

However, each of the summands in square brackets is by assumption 
of our claim as it applies to k — u < k. So any potentially non-zero terms 
contain factors of the form Y~\ with a ^ b only. But all such terms have 
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the same sign as (— l) fc - so that no cancellation is possible. Hence T^ k 
is either or has the sign of (— l) k . This is (I2.15p . 

Let E = {Eij} be the eigenvector matrix of T" 1 and U be its matrix 
of eigenvalues. 

To obtain (12.1 4p we claim that (T~ k )ij must be non-zero for some 
k < n; otherwise the eigenvectors of T will have to satisfy a homogeneous 
system of n + 1 linear equations: 

n 

^2 E i,kE kJ = 

k=l 
n 

Ei^Ekju^ 1 = 

k=l 

n 

^ E^Ekju^ 2 = 



fe=l 



J^^kE.ju^ = 0. (2.19) 

k=l 

Let E = (E i} iEij, Ei^E 2 j, ....E^ n E n j) and write (I2.19p as: 

UE T = 

where U is a square matrix with Uij = uj 1 ^ 1 for 1 < i,j ' < n. 

We know that E is a non-zero vector - otherwise Tij would be 0. 
Hence det U must be - that is it must have linearly dependent rows. 
Let 

u= fa 1 ,uZ\...u- x ) 

and 

u k = (u^ k ,u^ k ,...u- k ) 

Hence u fc_1 is the k-th row of U. The linear dependence of these vectors 
then implies that for some real numbers a^; for < k < n — 1, 

n-1 



fc=0 
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Which is, writing the vector form explicitly: For all i, 

akU- = 0. 

k=0 

So each of eigenvectors of T _1 is the solution to the same polynomial of 
degree n — 1. 

By the Fundamental Theorem of Algebra that would mean that at 
least two of our eigenvectors are identical. WLOG, let those be u\ and 
« 2 - Let E = (E itl E lt j +E i)2 E 2 j, ....E i>n E n j). Applying now the reasoning 
exactly identical to the above to the n — 1 dimensional vector E and 
matrix U, such that U^j = uj]^ 1 for 1 < i,j < n — 1 conclude that 
Uj = ui for some j ^ I. Proceeding by induction conclude that all 
eigenvalues of V are equal. This contradicts the irreducibility of V. 

Hence we know that ( 12. 1 5[) is indeed non-zero for some k < n. 

□ 

Proof of Lemma 11.21 concluded. 

Using this result to substitute into ( 12 . 1 8[) we conclude that whenever 
r-J = V/ < k and T~j 0, 

%j = ^r-j + o(t-^) = I|r-f | + o(r( fe+1 )), (2.20) 

and so is positive for t sufficiently large. 

□ 

We will say that has a zero of order h at if Vr < h 
rrj = and + 0. 

Remark 2.1. Note that the results of Lemmas \2j^ Theorem M . 31 and 
Corollary \2.1\ also hold for matrices that are not strictly positive definite. 

3 Upper bound on the critical point for 
Gaussian vectors. 

Proof of Theorem 11.41 and Corollary 11.11 To obtain the upper 
bound on ao from the Theorem 11.31 and Lemma 11.21 we use Lemma ll.lt 



12 



To apply the Lemma 11.11 we extract the coefficient of a special type of 
term in the following 



mi mj 



Lemma 3.1. Let A mijm . and B mijm . be the coefficients of the term s™*s 
(i 7^ j) in $i(S f ) and ^(S) respectively. Then for all rrii, rrij > 1 



j 



A„ 3 = l(^f)\l + 0(l/t)) (3.1) 



and 

„2 



= ^ Dl D 3 {-T^){l + 0(l/t)) (3.2) 
whenever T^j ^ 0. 

Proof Let m = mi + rrij. It is clear from (|1.5p that the only contribution 

to A miim ,. comes from 

(3.3) 



{(Qsr^ + nQsr}^ 



rn 

and, furthermore, we must have 



rrii rrij J_ J_ Qpo ,pi Qpi ,P2 %i-i,Pi "pm-i,Pm K^-^J 

Si S 'J 1=1 

with po = Pm = i', and all the other pi, 1 < I < m — 1, must be either i 
or j. Since 

rr 1 

g M = l + 0(l/t), g„ = l + 0(l/t) and %J = q hi = ^ + 0(l/t 2 ), 

(3.5) 

the terms on the right-hand side of fl3.4[) that are not 0(l/t 3 ) are those 
terms in which q it j and g^j each occur only once. This can happen in the 
following rrii ways 

Q r i,iQi,jQj,jQj,iQi^ r = 0, . . . , mi (3.6) 
(and, obviously, s = rrij — 1 and u = rrii — r) . Clearly 

vimvifliMti = %^ + WO- (3.7) 
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Repeating this argument with i and j interchanged we see that there 
are rrij ways that 



i(QS) n hj _ ( r J) 2 



0(l/t 3 ). (3.S 



Using (I3.4p - (I3.8I) we get (13.1 1) It is clear from (II. 6p that the only contri- 
bution to the leading term of B mum . comes from 



It 



Y,D u {q u MQS) m - 1 }i, j D j + Y,Du{qu, j s j {(QS) m - 1 } j , i D i 



,u=l u=l 



(3.9) 

By (13.51) we see that (I3.9P equals to 

^ (DiUQSrhjDj + DAiiQST}^) (1 + 0(l/t)); (3.10) 

where the contribution to the first summand of B mum . will be from a 
coefficient arising from m 8 of S; factors followed by rrij of s^- factors; and 
to the the second - from rrij of Sj factors followed by rrii of S{ factors. 
Furthermore, by (13. 5p the part of B mum . arising from {(QS) m }ij equals 
to 

flr^QB+OW*) (3.11) 
and similarly with i and j interchanged. Since 

CW^^^ (3 ' 12) 
we get (JO}. □ 

Proof of Theorem 11.41 continued. 

To begin, assume that T^ 1 ^ for all 
Let 



From (22) and (23): 



Cij{a) :=C mi , mj = [ -^ + " 2 Ai?i ! ( I +0(1//)) (3.13) 



t 1 
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Note that since the length of the term is fixed, for any N such that 
rrii,mj < N (13.131) holds for all t sufficiently large (i.e. t » N) thus 
making Cij(a) independent of rrii,mj. 

For Cij(a) to remain non- negative when DiDj < as t — > oo we need 
to have 

irr 1 ! 

^^ + a 2 D i D j >0. (3.14) 

Therefore, using Lemma 11.11 we get (11.81) . which is Theorem 11.41 for 
inverse covariance matrices with non-zero entries. 

Removing restriction T 1 ^ 

Now returning to (13. lft and ( 13.21) we see that if T _1 has a zero of 
order k at (i, j) for some < k < n — 1 then using Theorem 11.31 again 
we obtain: 

A ij = ^TT) (( r ~%) 2 + 0(r^) (3.15) 

And: 

B id = a 2 AjD ._i_| (r - fc) . .| + (t-^) (3.16) 

So that 

a» = ((r- fc )„) 2 + a 2 A/^|(r-%|(i + o(i/t)) (3.17) 



will be positive iff DiDj > 0, i.e. whenever z-th and j-th rows have 
the same sign. This concludes the proof of Theorem 11.41 and Corollary 
PI □ 



4 Elementary proof of Theorem 11.2 



We are now in the position to give the "'elementary'" - i.e. not involving 
the Markov process theory - proof of the Theorem ll.2[ This proof sim- 
plifies the proofs given in |4], Theorem 1.3; see also [3], Theorem 13.3.1. 
We will start by showing the result for the special case c = 1, see jl], 
Theorem 1.2: 
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Theorem 4.1. Let G be a mean zero Gaussian vector with strictly pos- 
itive definite covariance matrix T. The following are equivalent: 

1. (G + a) has infinitely divisible squares for all a G R. 

2. r _1 is an M-matrix with non-negative row sums. 
Proof 

(1 =>• 2): Let (G + a) have infinitely divisible squares for all a G R. 
By Theorem 13.2.1, [3j T _1 is an M-matrix. By Corollary D = in this 
case. Since T, and so also T _1 is positive definite this implies that all 
row sums of T _1 are non-negative. 

(2 1): Let T -1 be an M-matrix with non-negative row sums. Its 
Laplace Transform is given by Theorem [L3J The positivity of coefficients 
in $i follows from T -1 being an M-matrix; and that of coefficients in $2 
from Lemma 11.21 and row sums positivity. 

Proof of Theorem 11.21 

We first observe that if (1), (2) or (3) of Theorem 11.21 hold then 

• All components Cj of vector c must have the same sign. 

• If Cj = for some i, then c = 0. 

If (1) holds then no component of c is or c = 0. Suppose q = 
0, Cj 7^ or Ci,Cj have different signs for some pair Since 
(r + ac) 2 is infinitely divisible for all a, so is the sub vector (Tj + 
aci, Tj+acj). But for this sub vector r _1 Chas row sums of different 
sign. Hence it is not infinitely divisible for all a by Remark??. 
Hence Cj 7^ for all i and all components of c are of the same sign. 

If (2) holds then it also holds for any sub vector (Tj + q??, Tj + 
Cji], However by the same reasoning as above either r^ 3 or Tj^ 
is negative by the Lemma 2.1, [4]. 

If (3) holds then recalling that T _1 is strictly positive definite M- 
matrix conclude that q's are non-zero and of the same sign. 

The proof of Theorem 1.3 in [4] now shows that the result of 11.21 
follows immediately by replacing G by — . 
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5 Some other applications. 

Another open question pertaining to infinite divisibility of the Gaussian 
squares is: Suppose G 2 is not infinitely divisible. What could be said 
about infinite divisibility of (G + c) 2 ? By the technique employed in 
Theorem 11.21 it is enough to consider this question for (G + al) 2 for a e 
R. Furthermore, by conjugating G by an appropriate signature matrix 
one can consider vector G with covariance T having all non-negative 
entries. Since in two dimensions the inverse of covariance matrix is always 
an M-matrix, one should be concerned with vectors in three or more 
dimensions. 

Hence the question stands as follows: consider the gaussian vector G 
of dimension 3 or higher and assume that G 2 is not infinitely divisible. 
Let T be the covariance matrix of G such that Vij > for all(i, j). When 
is {G + a) 2 infinitely divisible for all al For no values of al Are there 
some vectors possessing the critical point - i.e. some a > such that 
(G + a) 2 is not infinitely divisible for a < «o and infinitely divisible for 
a > aol 

Theorem 11.31 allows us to shed some light on the above question for 
some special cases. 

Lemma 5.1. Let G be a Gaussian vector with positive, strictly positive 
definite covariance V. If there exist a triple i,j,k such that 



r 





(5.1) 
(5.2) 



and D i: Dj, D k > 



then 



(G + a) 2 is not infinitely divisible for any a G R. 



Proof 



Consider the coefficient of the term SiSj\ 




%j )(l + 0(l/t)) 



(5.3) 



This is negative when 




(5.4) 
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Now consider the term SiSjS k - Its coefficient is: 

a 2 

{Qi,jQj,k<lk,i + -j-[\ D i D Mj,k<lk,i + \DjD k \q itj q k>i + \D k Di\q id q^ k ]} 

•(1 + 0(1/*)) 

£3 l 1 1 j,k l k,i ^ 

^ ihn, vjy,;, - i\,ihY;]Y,:, - i),j),\\jyj:}{\ + o(i/t)) 

(5.5) 

Note that only the first term in square brackets is positive, while the 
other three terms are negative. Hence for this expression to be negative 
it is enough to have: 

« 2 < ||- < 5 - 6 > 

This shows that the series expansion of $ will have at least one neg- 
ative (non-constant term) coefficient for each a G R. 

□ 



Corollary 5.1. Let n = 3, 4 or 5. Assume that T > 0, is strictly positive 
definite and Y^ 1 has positive row sums. If Y^ 1 has at least one positive 
off-diagonal entry, then (G + a) 2 is not infinitely divisible for all a G R. 

In general case of Y, there exist a signature matrix iV such that Y = 
NTN > and replace r by Y throughout. 

Proof Note that to obtain this result it is enough to show that in 
dimension 3, 4 and 5 the inverse covariance matrix with positive off- 
diagonal elements has a triple of indices satisfying (15. ip . Also note that 
r _1 can not have a row with only positive entries. 

This is obvious for n = 3, since in this case one off diagonal element 
is positive, while two others are negative. 

To see the result for n = 4 and 5 assume that T _1 does not contain 
the required combination of entries. Note that in this case it is possible 
by means of row-column permutations to bring the matrix T _1 into the 
following form: 

if j > i and r^- > 0, then T~l > 0, for all k > j 
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When written in this form, it is easy to see that we must have T^ 1 = 
A + B, where A is the direct sum of M-matrices; B > whenever A = 0, 
and B = whenever A ^ 0. If n = 4 or 5 this is only possible if A is a 
direct sum of two such matrices, each of the size 2x2 for n = 4, or of size 
2x2 and 3x3 for n = 5. However, if this is the case, then we can find a 
signature matrix N such that iVT _1 iV is a weakly diagonally dominant 
matrix with negative off diagonal and positive diagonal elements. But 
the by [3], Remark 13.1.3, (iVT^iV) -1 = NTN has positive entries. 
Contradiction. Hence the triple of entries satisfying 15.11 indeed exists for 
n = 4 and 5. n 



6 General form of the coefficient for the 

term of the form s™ l s™ ] . 

1 j 

Theorem 1 1 . 3 1 can be used to derive the two index coefficient of $ explicitly 
for any dimension n. When applied to the two-dimensional case this 
result gives a shortcut for the derivation of expression for Rij, p used by 
Marcus and Rosen in [5]. 

In the statement of Lemma [37T1 we define A munij and B mum . to be the 
coefficients of the term s^s™ 3 ; {i ^ j) in ^(S 1 ) and Q^S) respectively. 
In (13. ip we estimate A mum ^ for fixed m^, rrij as t — > oo. However, to 
apply Lemma 11.11 to show coefficient positivity we must also consider 
Am^mj when m 8 and rrij are larger than t. We do this in the next lemma. 

Lemma 6.1. For rrii, rrij > 1 

min(mi,mj)-l , x . . 

5 ("7 XV) 



V + 



min(m,i,mj) — 1 2 

where 

'm; — l\ / m,- — 1 



nm,,™^) = (C l " (P+1) )(^" (P+1) )(^)( ; J ( % J- (6-2) 
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Proof By (II. 5p . to find A miTn . we sum over all terms in the trace 
of [QS) rrH+mi that involve Sj, m ; times and Sj, rrij times. Since the 
summands of the coefficient come from the trace of a matrix each of them 
will have the form qi 1 ,i 2 qi 2 ,i 3 ~~qi m _ 1 ,i m with i x = i m = i or i x = i m = j. 
It is convenient to consider the circular arrangement of the factors S(.), 
since once such arrangement is fixed each of m of its rotations result in 
the same coefficient. 

Consider arrangement of the s« and Sj into blocks. By blocks we 
mean an unbroken string of s< or Sj factors of length one or greater. For 
example, s^sfsj contains an Sj block of length 4 but SjS^sf does not. 

To account for all possible arrangements of indices of s^s™ 1 , consider 
any fixed circular arrangement. Suppose this arrangement has p + 1 
separate blocks of Sj, it then must also have p + 1 blocks of Sj. We will 
have < p < min{mj, rrij}. For each p there exist 

m ~ l \ h ~ l \ (6.3) 
p J \ p J 

ways to make this separation into blocks. Once the blocks are defined 
there exist exactly one circular arrangement of blocks( if we agree which 
index to always place first). 

Each arrangement can be rotated one position clockwise exactly m = 
mi + rrij times. A given arrangement will repeat itself after k one- 
positional rotations, where m = kv\ for some integer v%. Obviously 
also p + 1 = kv2 for some integer t> 2 . Hence this arrangement can be 
rotated v x times giving rise to a new arrangement. This cancels the V\ 
factor in the denominator of f ll.5p . leaving k. The coefficients in $x of all 
such rotations are identical and equal to 

To see this note that after any such rotation the resulting string of factors 
will either begin and end with the same index or will begin with Sj, end 
with Sj or vice versa. Using the symmetry of Q, see that: 

If it begins and ends with, say, Sj, the coefficient in trace (QS) m is: 

and is it begins with, say, Sj and ends withs,-, the coefficient in trace(Q5') m 
is: 

" V,)". 
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and the same with i,j interchanged. Using the fact that qij = qj^ again, 
see that this equals to the expression in figure parentheses above, although 
the order of q factors varies with the rotation; for example starting with 
an arrangement giving rise to a coefficient 

1 U\ VI v p+l 

rotating it one position results in a coefficient 

and so on. Also note that each fixed arrangement with a fixed starting 
point gives rise to two coefficient summands, each corresponding to the 
direction - clockwise or counterclockwise - of the order of g-factors. For 
example, the last arrangement '"read"' in the opposite direction will give 
the coefficient 

Vp+l Vl U\ 

m 

Therefore the coefficient of the whole sum is multiplied by 2. 

Finally we note that each of the v\ rotated block arrangements give 
rise to v 2 repetitions, according to which block is taken as the starting 
point of the circular arrangement. Therefore we divide by v 2 - Thus the 
resulting denominator is: v 2 k — p + 1. □ 

Remark 6.1. We pause in the development of the two-index coefficients 
to remark that the very argument just employed to find the denomina- 
tor of A mu „ Lj part of the coefficient is responsible for the non-existence 
of the nice '"closed" ' form for the coefficient in tree or more indices. 
Consider the coefficient of s™" 1 s • 3 s™ k in the power series expansion of 
$i(S f ), and assume that mi,mj,mk > 0. Let our factors be broken 
into Pi,Pj,Pk blocks respectively and assume an circular arrangement is 
fixed such that no two blocks of the same index are adjacent. We can 
rotate this arrangement u times to get a new arrangement, and obvi- 
ously m = TOj + rrij + m k = uv\ for some integer v±. Clearly, too, 
Pi = uvi,pj = uvj,pk = uvk- Now divide by m and multiply by u to get v\ 
in denominator as before. Now to account for different starting blocks we 
would have to divide by Vi,Vj or v k , depending on which block sequence 
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gives u new arrangements under rotation. But to know this we need addi- 
tional information on the blocks breakdown! This is what will prevent us 
from getting the sought coefficient based on the number of blocks alone. 

Recall that B m ^ mj is the coefficient of s^s™ 3 in the power series 
expansion of $2 (SO- Let 

mm(m,,mj)-l 

B mim = E /*„,,„.,, 

p=0 

Lemma 6.2. For irii,mj > 1 

min(mi,mj) — l ^ 

S mi , mj = —V(m i ,m j ,p){2D i D j 

Proof By (TOIL Q, (BHD and fT2TT2"j) we obtain the coefficient in P(£) 
arising from a particular arrangement, say s& 1 Sfc 2 s& 3 ...Sfc m . +m . , where fc^ 
is either % or j equals to: 



1 n 

- 2j ^r9(r,fci)?(fc 1 ,fc 2 )9(fe2,fc3) ' ' ' 
r=l 

+ 0(l/t)) 

Using see that (16. 5p equals to: 
-D kl q( kl ,k 2 )q(k 2 ,k 3 ) ■ ■ ■ 

To find B mum . jP , assume that one group of factors - either Sj or - 
is broken into p + 1 blocks for some p . Since we are dealing with a term 
in two indices only, one of the two cases is possible: 

1. Both groups of factors are broken into p + 1 blocks. Then < 
p < min{mj,mj}. There will be [ mi ~ )( mj ~ ) such arrangements. 
That means either starting with Sj, ending with Sj block, or vice 
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versa. Both will result in exactly 2p + 1 of q^j factors and have the 
coefficient in $ 2 ; 



min(m,i,m,j) — l 

E 

p=0 



rrii-1 
P 



rrij — 1\ a 2 



P J 2t 



A(CC) *'' J 



(2p+l) 



(2p+l) 



+^(co (g :^ )(P+1 ) ^Aj (i + o(vo) 



min(mi,m,j) — 1 g 

£ |-A^K,m„ P )(^ + ^)(l + 0(l/t)) (6.5) 

p=0 zr 9jj 



2. Now restrict p to 1 < p < min{mj, rrij}. The group of Sj factors 
is broken into p + 1 blocks and that of Sj factors into p blocks, or 
vice versa. (If rrii < rrij then there is also a case when there are m, 
of Si factors and rrii + 1 of Sj factors. The resulting coefficient 
in $ 2 is of the form B mumj , mi (l/t).) There will be p^" 1 ) ("i" 1 ) 

arrangements of the first kind and (^l^ 1 ) (""^'j of the second. This 
will result in exactly 2p of q^j factors, moreover the D factors will 
have the same index in each of the arrangements - same index as 
the factors broken into p + 1 blocks. The corresponding coefficient 
in $ 2 is, then: 



2 min(mi,mj)—l 

— Y (q mi q mj )q 2p 

2f L^i yHl > 1 H h3 > H h3 



p=l 



rrii — 1\ f m j ~ 1 



P 



P-I J (qi,iY p+1 KQjj) p 



+ 



rrij — 1 \ / m, ; — 1 



p — 1 



P /(?m) p (^) (p+1)J 



1 + 0(1/0) 



(6.6) 



Using now the identity: 



m — 1 
p — 1 



p fm — 1 



m — p \ p 
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see that (I6.6P equals to: 

min(m,i,mj) — l ^ 

%V{jn i , m h p) [-^D 2 + — - — Dj] (1 + 0(1/*)) 
z — ' It L m ; — p J rrij — p J 

(6.7) 

since the expression in square parentheses is when p = 0. 
Now denote as follows: 



q 2 

A mi>mj!P = V(m i ,m j ,p)-^j{l + 0{l/t)) (b\8) 



And 



,.2 



n "i.-/ , rr P D 2 ] P 

nii — p 3 TTLj - 

+2D i D J }(l + 0(l/t)) (6.9) 



B mi , mj:P = —Vim^m^p)^ D 2 + D 2 ] 

J It ^ L rrii—p J iJij — p J 



Let 



Hence, 



Crm,mj ,p -^mi,mj,p B mi m ^ p (6.10) 



min(m,i,mj)— 1 

Cmi,irij = ^ ] C mum:j) p (6.11) 

p=0 



And 



C. 



m,i,mj,p 



V(m,i,rnj,p) 



p+l 



a 2 
+ 2t 



2D i D j q iJ 



P 



TTLj — p 



TTli — p 



(6.12) 

;i + o(i/0) 



Denote the expression in figure parentheses by R mum ^ p . We are inter- 
ested in the sign of C mum ., which will coincide with the sign of C m . irn . iP , 
if they are identical for all p. Since the sign of Cmi,mj ,p is the same as 
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the sign of R mum )P we will concentrate on the later. Assume that m(mj 
is bounded by linear multiple of t 2 , i.e. 

^Jmiirtj < Nt for some fixed N 



R 



> 



> 



irr j 



" 2 q 2 r 2|A^-||r-/| 

, D2 i p 



p+1 t 2 '2 

A 2 p 



+ ( : 



t rrij — p t rrii — p' 



) }(i + o(i/0) 



ir 



-1 1 2 



'•J 



p+1 t 2 



+ 



a 



2i n.n.|ir _1 l 



2 L 

2\DiDj\ 



t 2 



P 



1 



ir; 



+ y|2A^|(-|rr/| + £)] (i + o(i/0) 



1 + 0(1/*)) 

p 



p + 



(6.13) 



The representation in ( 16 . 1 3f) easily leads to the proof of Lemmas 
5.2 and 5.3 of [5] by providing a combinatorial alternative to algebraic 
derivation of the expression for Rj,k,p employed in that paper. 



7 Some numerical examples. 

The Gaussian vector is called an assosiated with a Borel process when 
the covariance matrix is also a O-potential density of Borel right process. 
The following (necessary but not sufficient) property of an associated 
Gaussian vector is the direct consequence of the property of O-potential 
density of Borel right process: 

IV, <r,,Ar A/ . (7.i) 

contrasting with the regular property of the strictly positive definite co- 
variance matrix: 

it, < r M rv (7.2) 

Although the above regularity property is too weak to provide much 
information about G in general; in case n = 3 we have the following 
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Corollary 7.1. Let G = {G\,G2,G 3 ) be a Gaussian vector satisfying 
( [7. Then the upper bound for the critical point of G is non-zero. 

Proof This follows immediately from Theorem ll.4l by writing the matrix 
T in terms of T _1 and above mentioned regularity condition: 

Suppose r _1 has a zero entry 1^*- while Y satisfies (17. ip . Let 1 = 
gij and let 5-1,3 = 0. Then |ri j2 | = 01,203,3; |r 2 , 3 | = 01,102,3; ^2,2! = 
01,103,3- This implies that g 1<2 < gi,i and #2,3 < 93,3- Hence D x and D 3 
are both positive, which makes the first case of Corollary H] impossible 
for such matrices. 

The G- vectors with inverse covariance matrices T _1 with zero entries 
and corresponding row sums of different signs indeed exist. 
Below is an example for n = 3 with one zero entry: 

Let 

1 / 15 4 2 \ / 1 -2 \ 

T = - 4 2 1 ; then T^ 1 = -2 8 -1 
7 \214/ \ -1 2 / 

Where T -1 is an M-matrix, so G 2 is infinitely divisible and T^ 2 = 2 
and T^j = 0, D 1 = -1 < 0, D 3 = 1 > 0. 

Therefore for a vector with this covariance matrix = —2^- while 
Aij = h so that Cij = 4 — 2a 2 i is negative for all sufficiently large t 
and all a 7^ 0, so the corresponding vector has critical point a = 0. Note 
that T does not satisfy (17. ip . 

Needless to say, examples of such matrices abound in higher dimen- 
sions. 

The following example illustrates that (17. ip is not a sufficient condi- 
tion for a vector to be associated: 




Where T satisfies (17. ip . T -1 is an M-matrix and has the first row sum 
negative - hence our vector is is not associated. 

It is possible, of course, for a 3-dimensional vector to satisfy (17. ip and 
have all positive row sums - i.e. be an associated vector and still have a 
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zero entry in the inverse covariance matrix: 



. / 15 4 2 \ / 3 -2 \ 

r = — 4 6 3 ; then T" 1 = -2 8 -1 
6( \ 2 3 20 / \ -1 2 / 

In dimensions 4 and higher the condition (17. ip no longer prevents the 
r _1 from having a zero entry at the intersection on the row sums with 
different signs: 



( 10 


-3 


-3 


9 


-3 


-2 


I o 


-3 



-3 \ 

-2 -3 

9 -3 
-3 



65 
11 



, then T 



257 


39 


39 


99 


1400 
39 


280 
171 


280 
115 


700 

33 


280 
39 


616 
115 


616 
171 


140 

33 


280 
on 


616 
33 


616 
33 


140 
1 13 



700 



140 140 



350 



and satisfies (17.11) . 

An example below is of the 3-dimensional Gaussian vector with pos- 
itive, strictly positive definite covariance matrix T and the diagonally 
dominant with one positive off diagonal entry. Hence the set of 
Gaussian vectors in Lemma 15.11 is non- vacuous already in dimension 3. 

1 / 15 4 2 \ / 23 -14 4 \ 

r = — 4 8 5; then = -14 86 -67 

297 \ 2 5 6 / \ 4 -67 104 / 
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